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Abstract 

The randomly pinned planar flux line array is supposed to show a phase 
transition to a vortex glass phase at low temperatures. This transition has 
been examined by using a mapping onto a 2D XY-model with random aniso- 
tropy but without vortices and applying a renormalization group treatment 
to the replicated Hamiltonian based on the mapping to a Coulomb gas of 
vector charges. This renormalization group approach is extended by deriving 
renormalization group flow equations which take into account the possibility 
of a one-step replica symmetry breaking. It is shown that the renormalization 
group flow is unstable with respect to replica asymmetric perturbations and 
new fixed points with a broken replica symmetry are obtained. Approaching 
these fixed points the system can optimize its free energy contributions from 
fluctuations on large length scales; an optimal block size parameter m can be 
found. Correlation functions for the case of a broken replica symmetry can be 
calculated. We obtain both correlations diverging as In r and In 2 r depending 
on the choice of m. 



1 Introduction 

The technological aspects of high-T c superconductors in strong magnetic fields and 
especially of their ability to preserve superconductivity by flux pinning |IJ have 
led to intense theoretical studies of the properties of a flux line array in a type-II 
superconductor with random point-like pinning centers [Q, ||, ||, Q [L(], [12], [14], |I~6|| . 



It has been conjectured 0, ||, £| that the flux lines in a superconductor with point 
disorder form a new thermodynamic phase, the vortex glass phase. It is supposed 
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that in this phase the flux lines are collectively pinned by the point defects and 
energy barriers between different metastable states of the flux line array occur which 
diverge with increasing length scale L leading to a glassy dynamics and zero linear 
resistivity jl2|, [15], [16|] . But there is still too little conclusive evidence to confirm this 
scenario by analytic means. 

Due to the absence of topological defects [15], [Kj in 1 + 1 dimensions, the planar 
1 + 1-dimensional flux line array can be well treated analytically in an elastic ap- 
proach. Actually, the system of flux lines in a type-II superconducting plane with 
parallel magnetic field and point disorder is the only system for which the existence 
of a vortex glass phase has been shown analytically 0, [51 0, |10[ |12[ by applying 
various different methods. 

On the other hand, the predictions for important physical features of the disor- 
dered 1 + 1 dimensional flux line array obtained by the different analytical methods 
differ significantly. In particular, there are still many competing conjectures con- 
cerning the correlations in the vortex glass phase. Moreover, the results of numerical 
simulations confirm neither of the analytical predictions |14|]- Essentially three ana- 
lytic approaches have been applied to the problem: (i) After using the replica trick 
and mapping onto a 2D XY-model with random anisotropy but without vortices 
[0, 0, a renormalization group (RG) calculation || 0] has been carried out with 
the replicated Hamiltonian not taking into account the possibility of replica sym- 
metry breaking (RSB). (ii) The replicated Hamiltonian has also been studied by a 
variational treatment admitting of continuous RSB finding that a one-step breaking 
is realized ||, [H|. (iii) Without making use of the replica trick, the corresponding 



kinetic equation has been treated by a dynamical RG analysis [11, 12 



In the present paper we want to study how the concept of RSB could enter into 
a RG analysis. For this purpose we map the disordered planar flux line array onto 
the 2D XY-model with random anisotropy and perform a RG calculation with the 
replicated Hamiltonian where we generalize the set of coupling constants such that 
we can take into account a one-step RSB. Our aim is to show that an instability 
with respect to one-step RSB can also be found in the RG treatment. This leads 
to a more unified view of the approaches (i) and (ii), and results obtained by the 
variational approach can partly be reproduced in our calculation. 

On scales larger than the flux line distance / the planar flux line array is described 
by an elastic model with the positions of the flux lines given by a scalar transverse 
displacement field 0(r)//7r. The field 0(r) itself can be regarded as a phase field 
giving the phase shift of the superconducting phase caused by the flux line displace- 
ments, because an increase in the flux line displacement by I induces a shift of ir in 
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the field <p(r). On large length scales, the planar flux line array interacting with the 
random pinning centers can then be described by the Hamiltonian [0, || [5], [17]] 

Ln[<f>] = J d 2 r {^(V0) 2 + Vi{r) sin(20(r)) + V 2 (r) cos(20(r))| . (1) 

The second term containing the random potential V(r) with zero average and short 
range correlations (i,j = 1,2) 



^(r)^(r0 = 2^(r -r')<% (2) 

models the interaction of the flux lines with the point disorder in the continuum 
description (g includes a factor 1/T 2 ). The function 5^ is a delta-like function of the 
small width £ given by the maximum of the flux line core radius and the impurity 
size. A crucial feature of this term is to respect the periodicity of the lattice, i.e. it 
is invariant under a uniform shift ir of the field 0(r). By rescaling of one coordinate 
the isotropic first term for the elastic energy is obtained with one elastic constant 
K (including a factor 1/T). 

Each of the three approaches sketched above 0, [5], ||, [7j |9|, [10], [11], [12| as well as 



the numerical simulation fl4[] , yield a phase transition ati^ = l/7rorr = with 



which serves as a small parameter measuring the distance from the transition and 
controlling expansions around the transition. For r < the system is in a high- 
temperature phase, disorder is not relevant and does not alter the correlations in- 
duced by thermal fluctuations on large length scales (((f>(r) — (f)(0)) 2 ) = (lnr) / '(ttK). 

However, the results concerning the correlations for r > in the glassy phase 
differ significantly. The RG analysis carried out on the replica Hamiltonian in a 
replica symmetric way yields correlations ((<f)(r) — 0(O)) 2 ) oc r 2 ln 2 r at large length 
scales r [7J, §). The same result is obtained in the dynamical RG calculation 
[p] , [T2|l . Besides, correlations diverging like the square of the logarithm follow from 



a real-space RG procedure [I3"|. On the other hand, in the variational approach with 



one-step RSB, correlations ((0(r) — (f)(0)) 2 } = (lnr)/(7r_ft'(l — r))) are found to 
diverge logarithmically but with a prefactor increasing with decreasing temperature. 
Logarithmically diverging correlations have also been found in the numerical sim- 
ulation but the prefactor of the logarithm does not accord with the analytical 
prediction in M. Our calculation including one-step RSB in the RG analysis of 



3 



the replicated Hamiltonian can reproduce both correlations diverging with a simple 
logarithm and correlations diverging with the square of the logarithm depending on 
the choice of the block size parameter m in the one-step RSB scheme. 



2 RG Analysis 

Introducing n replicas and averaging over the disorder gives the effective replicated 
Hamiltonian (with replica indices a,/3 running from 1 to n) 

hi R [<t>*\ = Jd*r\W K^V^ • - ]T 9ap cos (2(0 Q - <j> p )) \ (4) 

\ a/3 a/3 ) 

with matrices K a p and g a p taking on their bare values 

K a p,o = K5 a/ 3 (5) 
9a0fl = 9 ■ (6) 

This Hamiltonian is equivalent to the replica Hamiltonian of a 2D XY-model with 
random anisotropy but without vortices ||. The variational studies allowing for 
continuous RSB performed so far on this Hamiltonian |], [H| have shown that in the 
2-dimensional model considered here a one-step RSB is realized. Our calculation 
is restricted so far to a one-step RSB scheme but the results from the variational 
approach suggest that our results may stay valid even if it is possible to extend the 
calculation to higher steps of RSB or a continuous RSB scheme. 

Introducing one-step RSB and following the resulting RG flow, it is necessary 
to admit matrices K a/3 and g a p of the form K a p = A8 a p + B5 a p + C and g a/ 3 = 
9i8 a p + #2(1 — (W); the elements of the matrix 5 a/3 are 1 if a and f3 belong to the 
same block of size m and otherwise. Because of (0), (|J) we have initially A = K, 
Bo = C = and gi = #2,0 = 9- To derivate the full RG flow equations, we perform 
an analysis technically similar to that of Cardy and Ostlund [BJ but with significant 
extensions to take into account the one-step RSB. The calculation is based on the 
mapping onto a coupled Coulomb gas. We want to consider only weak disorder, so 
initially the disorder strength g will be small; also throughout the RG procedure 
the matrix elements of g a p stay sufficiently small to use standard methods || |H| 
to transform the cos-couplings in the partition sum and to integrate out the fields 
4> a (r) in favour of integer charges n^tr) (a < (3) with fugacity g a p. The replicated 
disorder averaged partition sum Z n factors then into Z n = Z e \Zc where the factor 
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Z e i represents the purely elastic part of the partition sum; this factor plays a role only 
in deriving the RG equation for the free energy density and will be considered later 
on in detail. In the Coulomb gas factor Zc of the partition sum, it has to be summed 
over all spatial configurations of interacting charges, which can take on any integer 
value, but because the charge fugacities g a @ are sufficiently small, only positive and 
negative unit charges have to be considered, which obey in addition a neutrality 
condition. Switching from the continuum description with a short wavelength cutoff 
£ to a description on a square lattice with lattice constant £ and lattice vectors R 
for easier notation, one obtains the following expression for the replicated disorder 
averaged partition sum: 

= Z el Z c = Z el x\{\{ E ex P ("Wo) (7) 

R a</3 n Q(3 (R)=-l 

-He = \ E E W «/?( R ) M R ') G'(R-R') + 

* R^R' a,/3,7,<5 

+ EE MR)) 2 lng aP , (8) 

R a/3 

where G'(R) = ln(|R|/£) and np a := — n a p (a < (3). The matrix (K^ 1 )^ = 
a5 a /3 + bS a /3 + c has the same block form as K a p with a = l/A, b = —B/A(A + mB) 
in the limit n —>■ 0. 

The block form of the matrix g a p implies that two kinds of charges exist differing 
in their fugacities and, moreover, in their interactions with other charges due to the 
block form of the matrix (K^ 1 )^. 

Taking this into account, a RG calculation in the style of Cardy and Ostlund 
(CO) |0 can be performed, which yields the following RG recursions in the limit 
n — > upon a change of scale by a factor e l (Henceforth we always include a factor 
2/V in a, b and c and a factor 4-7r£ 2 in g 1; g 2 .): 



da/dl = — | a 2 m (g\ — gl) (9) 

db/dl = | a 2 (gl-g 2 2 ) (10) 

dc/dl = | Oq g\ (11) 

d 9l /dl = {2-a)g x - ±(2 - m)g\ - \mg\ (12) 

dg 2 /dl = (2 -a -b)g 2 - mgl - (l-m)g 1 g 2 . (13) 



The parameter m is a free parameter in these equations with < m < 1 in the limit 
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n — > 0; possible choices for m will be discussed later on. (|9]), ( [TQ|) show that 

a + m& = a = 2 - 2r (14) 

is not renormalized; this result is exact to all orders in the gi (i = 1, 2) due to a 
statistical invariance under tilt || ^j. 

In the special replica symmetric cases m = 1 and m = 0, we get back the flow 
equations of CO: In the RG equations for m — 1 (m = 0), g 2 (gx) plays the role of 
the single disorder strength parameter in ||, the diagonal matrix elements a + b (a) 
of K~p are not renormalized, and also the off-diagonal matrix elements c (b + c) and 
the fugacity g 2 (gi) renormalize as in ||. The system exhibits the known CO fixed 
points g 2 = (g{ = 0) and g 2 = 2r (g{ = 2r). For m = 1 the RG flow is sketched 
in Fig. la. g\ (g 2 ) does not feed back into the RG flow of the other quantities 
and does therefore not enter physical results like correlation functions (see below). 
For this reason the introduction of a small initial replica asymmetric perturbation 
Ag = <7i,q — #2,0 has no effect on physical results if m = 1 (m = 0), although Ag 
turns out to be a relevant perturbation under RG (see below). 

Starting as in (§) with replica symmetric initial conditions gx^ = ^,0; replica 
symmetry is preserved throughout the RG procedure independently of m, and a, b 
are not renormalized; therefore the CO scenario with the trivial fixed point g\ = 
g\ — and the non-trivial CO fixed point g\ — g\ — 2r is reproduced if replica 
symmetry holds initially. 

However, introducing a small initial replica asymmetry Ag = gi t o — (72,0 7^ 
contrary to (||), the RG flow develops for r > an instability with respect to RSB. 
The system flows for Ag > to a regime with g\ > g 2 and for Ag < to a regime 
9i < 92 (entering on large length scales the unphysical regime of negative fugacities 
gi). In particular the replica symmetric CO fixed point g\ = g\ = 2r, a* = ao is 
unstable against small replica asymmetric perturbations. A linear stability analysis 
of the CO fixed point yields (Aa = — a) 

dAa/dl = -ma^T Ag (15) 

dAg/dl = 2—t Aa . (16) 
m 

These equations describe an instability with eigenvalue 2(1 — r)r of the CO fixed 
point with respect to perturbations Ag. To avoid entering the unphysical regime 
of negative fugacities, we consider only perturbations Ag > 0. As it is seen from 
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(|T5D, (|T6|), such a perturbation causes the charge interaction strength parameter a 
to decrease and the asymmetry Ag to increase; finally, a renormalizes to following 
(|9|). This flow towards the fixed point a* = implies that one non- interacting type 
of unit charge with fugacity g\ appear on large length scales. Furthermore, we can 
find from (|l"2"D , (|13D two additional non-trivial RSB fixed points (|17D and (0) with 
a* = 0, b* = a /m for each of them: 



9i 



1 — m , 2 4 x1 /9 

2 a + a 1 + — 1/2 

m m a 



2 — m 1 — m . 2 4. 1/9 

g% = 2 a do (1 + — ) 1/2 and (17) 

m m m ao 

9a = , 9l = . (18) 

At a certain 

m * = 1 - r /3 + C(r 2 ) (19) 

the fixed points (0) and (|T8|) fall exactly together. Only for m* < m < 1 the fixed 
point (|17D is in the physical regime g\ > of non-negative fugacities. Moreover, the 
fixed point (p!7[) is in this range of m stable with respect to perturbations in g\ and #2 
(getting marginal with respect to perturbations in g2 at m = m* where it coincides 
with (|T8|)), whereas the fixed point ( pT8[ ) is unstable with respect to perturbations 



#2 > 0. Therefore the fixed point (|17D is attractive for all RG trajectories starting 



with gi o > (72.0 > ( as illustrated in Fig. lb) while the fixed point (IB) is attractive 



only for RG trajectories with #i j0 > g 2 ,o = 0. For < m < m* (0) is the only RSB 
fixed point in the physical regime of non-negative fugacities gi > 0. It is in this 
range of m the attractive fixed point for all RG trajectories with gi t0 > g 2 ,o > (see 
Fig. lc,ld); furthermore, it is stable with respect to perturbations in gi and g 2 . 

As pointed out in (|5[), (0), the proper initial values K a p tQ = K5 a p and g a /3,o = 9 
are replica symmetric with Ag = 0. It remains unclear in this approach how the 
initial asymmetry Ag > necessary for the development of an instability with 
respect to RSB can be obtained from physical reasons. One hint is given in the 
next section where it is shown by comparison with the replica symmetric CO flow 
for Ag = that contributions to the free energy from large scale fluctuations can 
be optimized (which means maximized in the limit n — > 0) if a small perturbation 
Ag > is introduced. 

In the high-temperature phase for r < the system flows to the stable trivial 
replica symmetric fixed point g* = g\ = regardless of an initial asymmetry Ag ^ 
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0. In this phase the trivial replica fixed point is stable with respect to the RSB 
perturbation Ag so that RSB cannot occur in the high-temperature phase as it is 
expected. For r = the trivial fixed point stays marginally stable. 

3 Free Energy and RSB 

We want to proceed with a discussion of energetic aspects of the instability in the 
RG flow upon introducing a replica asymmetric perturbation Ag > in the low- 
temperature phase. This enables us to fix the so far undetermined block size param- 
eter m if Ag > and to compare the free energy in the RSB case with Ag > with 
the free energy in the replica symmetric CO case Ag = 0. The standard procedure 
to determine m is to maximize the free energy density per replica in the limit n — > 
with respect to the additional free parameter m. 

In the RG approach it is possible to derive the RG flow equation for the free 
energy density and calculate the free energy by integrating the flow equation; more- 
over, one can separate the contributions to the free energy from fluctuations on 
different length scales examining the flow of the free energy. In the RG procedure 
of increasing the cutoff £ to £e dl in the partition sum and rescaling the scale, one 
collects contributions not renormalizing the coupling constants; these contributions 
enter the renormalization of the free energy. Such terms are generated both in 
the factor Zc by contributions from integrating out the charge configurations as 
described above and in the factor 

Zei = exp(-i(L/2vr0 2 /" /^ 2 g{^ln(£V/2vr) + 

+ ln ( 2 fl ) " - ln ( 2 ao) + T mo] j J (20) 

(in the limit n — > and with L denoting the linear dimension of the system) by 
contributions from increasing the cutoff £ and adjusting the couplings a, c according 
to the flow equations (|S|), flTT|). Finally, one obtains in the limit n — > for the free 
energy density per replica (apart from an additive constant independent of m) the 
RG recursion relation 

df/dl = 2f + —dZ/dl 

Ttlj 

= V ~ -^- 2 (mgl + (l-m)gl) - 



S 



From this recursion relation the initial free energy density /o can be obtained by 
following the flow: 



Using ([H]), (^) we want to study the contributions to the free energy from fluc- 
tuations on different length scales in the low-temperature phase for the flow to the 
RSB fixed points fllTD , (p~8|) induced by a small initial perturbation Ag Q > and, 
in comparison, for the replica symmetric flow to the CO fixed point starting with 
Ag = 0. To keep calculations tractable, we choose initial conditions in the vicinity 
of the CO fixed point, i.e. = 2r + Ag , #2,0 — 2r, when examining the RSB case. 
For general initial conditions with Ag > there is initially a flow towards the CO 
fixed point slightly perturbed by the small asymmetry Ag > 0; for the study of the 
energetic effects of the RSB instability in the flow, this essentially replica symmetric 
part of the RG flow should be negligible. 

As a consequence of the factor e~ 21 appearing in the integrand of (p2|), the main 
contribution to f comes from the short scales. Starting at g i>0 = 2r + Ag , g 2 ,o = 2r 
and evaluating fl2"2| ) straightforwardly to the leading order in Ag , one obtains a 
maximum of /o a t the replica symmetric m — 0. This is because the main con- 
tribution to —(dZ/dl)/nL 2 comes on short scales from the term [— 7r 2 a c/2]/2£ 2 . 
The m-dependent part of [— 7r 2 aoc/2]/2£ 2 can be approximated by means of a linear 
stability analysis of the flow equations ©-(El) & t t ne CO fixed point enlarging on 
([ID, (|T§ as [A# (7r 2 ajj/16)(l - m)(exp (2rZ) - l)]/2£ 2 with a maximum at m = 0. 

On the other hand, the maximization of the asymptotic large scale contributions 
leads to a quite different result. Examining the large scale contributions, one has to 
investigate the asymptotics of the integrand in (^) and to maximize — (dZ/dl)/nL 2 
in the limit of I —>■ 00. For this purpose it is necessary to derive the asymptotics of 
c which is determined by the flow equation (|TT|). From (|TTD follows that for m < 
m* < 1 c is asymptotically linear divergent with an asymptotics c{l) ~ / g^il— 1~) 2 /2 
where g% is taken in the stable RSB fixed point (|T7|), which is g^ — 6(m — m*) to 
a good approximation. However, we find in the regime < m < m* from (|TI|) a 
saturation of c to a value c* because the stable RSB fixed point is in this regime 
given by ( |18D with g% = 0. To obtain an estimate for c* one has to determine 




(22) 
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the characteristic scale I* on which gi renormalizes towards 0; a linear stability 
analysis of the flow equations (|T2"D, (p~3|) at the CO fixed point extending ([15]), 
QISD reveals that Z* can be approximated as I* ~ (ln(4r/(l — m)Ag ))/2r. From 
( |Pi~D it follows c* ~ rln(4r/(l — m)Ag Q ) for the leading order contribution in r. 
Moreover, it is seen from the flow equation @ that a(l) ~ 1/Z on large scales. Using 
these results for c and a, one can verify easily from (|21|) that the most divergent 
contributions to ~(dZ/dl)/nL 2 come from [— 2 In (cto/a)(l —m)/m — 7r 2 aoc/2]/2£ 2 
for large Z. Maximization of these terms yields 

m = m* = 1 - r/3 + C(r 2 ) , (23) 

because the maximization of the second term restricts m to values < m < m* to 
avoid the occurrence of the linear divergence in the regime m* < m < 1 and the 
maximization of the only logarithmically diverging first term singles out the greatest 
value m = m* of the interval < m < m*. (|23"D is in fairly good agreement with ||. 

RSB is a large scale effect associated with the existence of diverging energy 
barriers generating metastable states. Therefore it seems to be more reasonable to 
consider only the large scale contributions to the free energy in (j2~2"D although the 
expression Q2"2]) for f is dominated by its short scale part. This is equivalent to 
considering the free energy of the renormalized but not rescaled Hamiltonian on 
large scales but discarding a constant energy shift depending on m which comes 
from short scales. This energy shift, which is essentially replica symmetric, may 
describe the free energy of single metastable states. In the presence of an initial 
asymmetric perturbation Ag > 0, maximization of the large scale contributions to 
the free energy yields then a maximum at m = m* as derivated above. 

Comparison of these large scale contributions for the flow to the RSB fixed point 
when ( Ag > 0) and for the replica symmetric flow to the CO fixed point ( Ag = 0) 
shows that this part of the free energy is greater in the RSB case. This is because 
the most divergent contribution to —(dZ/dl)/nL 2 is in the replica symmetric case 
as in the RSB case given by [-7r 2 a c/2]/2£ 2 with c(l) ~ I gf(l - r) 2 /2 but the 
replica symmetric CO fixed point value 2r for g\ is always greater than or equal to 
(m = 1) the RSB fixed point values given by QT7D, flTBp. Therefore it is energetically 
favorable on large scales to break the replica symmetry by introducing a perturbation 
Ag Q > 0. This energy gain can occur on scales larger than L* = £ exp (l*\ m=m *) = 
£(4r/(l - m*)Ag )V*r = m/ Ag Y^. 

In the high-temperature phase the trivial fixed point is stable with respect to the 
introduction of a perturbation Ag ^ 0. For this reason the large scale contributions 
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to the free energy are the same as in the replica symmetric case. For the short 
scale contributions to the free energy the same argumentation applies as in the 
low-temperature phase leading to a maximum at the replica symmetric m = 1 if 
Ag 7^ 0. So RSB is energetically not favorable in the high-temperature phase as it 
is expected. 

4 Correlations 

The RG flow and fixed point structure changes significantly upon introducing an 
energetically favorable, replica asymmetric perturbation Ag > in the low-tem- 
perature phase as outlined above. As well the behaviour of the (00)-correlations 
changes drastically depending on the value of m. The Fourier transformed correla- 
tions between replicas on large scales can be calculated by using a Gaussian approx- 
imation to the renormalized but not rescaled replica Hamiltonian, which yields for 
small q ~ e~ l /£ g 

(0 Q (q)<M-q)> = j 2 (K- l U{l = \n^/q)) , (24) 

so that the large scale correlations depend essentially only on the asymptotic RG 
flow of the matrix elements a, b and c. 

From expression ( p4]) one can verify that the connected correlation function 

(0( q )0(- q ))-<0( q ))(0(-q)) = liml/n^^q^-q)) = (1 -r)n/q 2 (25) 

does not change its form at the transition, independently from the introduction of 
a nonzero Ag due to the non-renormalization of a + mb = a , contrary to the 
(00)-correlations 0. 

For the (00) -correlations ( p4|) yields 

(0(q)0(-q)> = (0 Q (q)0 Q (-q)> = ^ WO + 6(0 + c(Q)|wn«/ fl ) . (26) 

In the high-temperature phase no RSB takes place, even if Ag > 0, and there is 
essentially no renormalization of a, b and c ; the connected (00) — (0) (0)-correlation 
function and the (00)-correlation function coincide then and ((0(r) — 0(O)) 2 ) = 
(l-r)ln(r/0- 
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In the low-temperature phase the asymptotics of c, which is determined by the 
flow equation fllTD , is of special interest because in the replica symmetric case, i.e. 
without a replica asymmetric perturbation (Ag Q = 0), the asymptotics c ~ 2(1 — 
r) 2 r 2 / diverging linearly gives correlations ((0(r) — 0(O)) 2 ) ~ ((1— r) 2 r 2 /2) In 2 (r/£) 
diverging with ln 2 r ||. 

In the RSB case with an initial Ag > 0, c has also a linear divergent asymptotics 
c(l) ~ I <?2 2 (1 — r ) 2 /2 for m * < m < 1 (see above), where g% — 6(m — m*) is taken in 
the stable RSB fixed point fllTT). This entails (00)-correlations diverging also with 
In 2 r for this range of m but with a pref actor reduced by a factor 9(m — m*) 2 /r 2 < 1 
compared to the replica symmetric case; in particular, we get back the replica sym- 
metric CO result choosing m = 1. The situation changes significantly in the regime 
< m < m* because the stable RSB fixed point is in this regime given by ( |T8D with 
g>2 = 0. Therefore c saturates on large scales to a value c* ~ r In (4t/ (1 — m)Ag ) 
for the leading order contribution in r (see above) and we obtain from (^) only 
logarithmically divergent (</>0)-correlations with a prefactor a /2m + c*/2 which is 
greater than in the high-temperature phase. To the leading order in r this yields 
correlations ((0(r) - 0(O)) 2 ) ~ ((1 - r)/m + r (ln (4r/(l - m)A g Q ))/2) In (r/Q; 
with our above choice (|23D of m ~ 1 — r/3 we get ((0(r) — 0(O)) 2 ) ~ ((1 — 2r/3) + 
r(ln (12/A(7 ))/2) In (r/^). This implies that the prefactor of the logarithm increases 
with r in the low-temperature phase. 

Our results for the low-temperature phase show that within the one-step RSB 
RG approach with a small initial replica asymmetry Ag > 0, it is possible to 
obtain the known replica symmetric result for the (00)-correlations diverging like 
In r H if m = 1 as well as (00)-correlations diverging in the same way but with 
a smaller prefactor if m* < m < 1 and logarithmically divergent (0</>)-correlations 
with a prefactor increasing with decreasing temperature if < m < m*. The latter 



possibilities are of interest with regard to the numerical simulations |T4j] and the 
results of the variational approaches || |i~0] . 



5 Conclusion 

To summarize we have shown an instability in the RG flow of the disordered planar 
flux line array, which is equivalent to the 2D XY-model with random anisotropy 
but without vortices, with respect to a one-step RSB. The flow approaches new 
RSB fixed points if an initial replica asymmetric perturbation is introduced. The 
system can optimize its free energy contributions from large length scale (> L*) 
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fluctuations by breaking the replica symmetry and approaching the RSB fixed point 
where the energetical optimal choice of m is m = m* ~ 1 — r/3. Introducing the 
initial perturbation, the (00)-correlations show a In r-divergence on large length 
scales in the range m* < m < 1 returning to the replica symmetric result at m — 1; 
for < m < m* the correlations diverge only as lnr, which is especially for m = m* 
the case. 

During completion of this work P. Le Doussal and T. Giamarchi have submitted 
a letter |19[ in which they find independently from our results an instability with 



respect to RSB in the 2D XY model in a random field. 
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Figures Caption 

Fig. 1: RG flow trajectories for r = 0.05 and different values of m with an initial 
replica asymmetry Ag = (7i,o — ^2,0 — t2 ] a) m = 1, b) m = (l+m*)/2, c) m — m*, 
d) m = 1 — 5r. The dashed line is the line g\ = gi of replica symmetric values. 
The trivial fixed point and the CO fixed point on this line as well as the RSB fixed 
points flTTD , dl8|) are plotted. 
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